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Abstract 

A review is given on the recently proposed two dimensional axion model 
(0(3) fj- model with a dynamical 0-term) and the T-duality relating it 
to the SU(2)xU(l) symmetric anisotropic cr-model. Strong evidence is 
presented for the correctness of the proposed S'-matrix for both models 
comparing perturbative and Thermodynamical Bethe Ansatz calculations 
for different types of free energies. This also provides a very stringent 
test of the validity of T-duality transformation at the quantum level. The 
quantum non-integrability of the 0(3) cr-model with a constant ^-term, 
in contradistinction to the axion model, is illustrated by calculating the 
2 — > 3 particle production amplitude to lowest order in 9. 
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1 The axion model 



Let us consider the following two dimensional cr-model defined by the LagrangianQ: 

£ = Xdau'^d^'n^ + rd.ed^e + — e'^"e"''%'^an^a,n'= , (1.1) 

2A ^ 327r2(l + ^) ^ Stt m ^ > v ; 

where A, g are (real) parameters (couplings) and the n-fields are subject to the 
constraint n"n" = 1. I shall refer to ( |1 . 1| ) as the 'axion model' since it is the 
0(3) nonlinear cr-model with a dynamical ^-term, which is in turn a natural 
two-dimensional analogue of its phenomenologically important four- dimensional 
counterpart, in the same way as the 0(3) a-model is considered to be a 2D 
analogue of 4D non-abelian gauge theories. As will be exhibited later (|1.1| ) is 
asymptotically free for the coupling range — 1 < ^ < 0. The theory defined by 
(|1.1|) is an 0(3) non-linear cx-model, coupled to a scalar field, 6, (whose normal- 
ization has been chosen for later convenience) through the Hopf term. The latter 
is proportional to the topological current of the 0(3) model and with the normal- 
ization chosen in ( |1 . 1\ ) its space-time integral (after a Wick rotation) yields the 
topological charge, which is integer valued (for non singular configurations). The 
variable 6 therefore has a natural interpretation as an angle, 6 G [0,27r], so that 
the three dimensional target space of the axion model is topologically S"^ x S^. 
Although at the classical level only one of the parameters (couplings), (A, g), is 
relevant I prefer to introduce them from the very beginning as in the quantum 
theory both couplings play a role. It is not difficult to see that for ^ ^ — 1 the 6 
field decouples from n", i.e. one obtains the 0(3) a-model and a decoupled free 
scalar field. There is another special value of the coupling, ^ = 0, when ( |1 . 1| ) 
turns out to be canonically equivalent to the 0(4) a-model. 

The 'manifest' global symmetries of the axion model are S0(3) xSO(2)0, where 
SO(2)6) generates a shift 9^9 + const., a symmetry of (|Ll|) up to a total 
derivative. 

In this paper I shall review the axion model and its dual, the SU(2)xU(l) 
symmetric anisotropic principal cr-model. The most remarkable feature of the 
axion model is its conjectured quantum integrability (in the sense that there is 
no particle production). It allows one to deduce by bootstrap methods its exact 
spectrum, which turns out to depend on the value of the (renormalization group 
invariant) coupling ratio p = 27r(l + g)/X. As long as < p < 1 the spectrum 
consists of S0(3) triplets and singlets ('breathers') as well as of doublets ('kinks'), 
while for p > 1 only two doublets. 

This review is organized in the following way: In Section 2 it is shown that 
the classical equations of motion of the axion model admit a (standard) Lax pair, 

^The following conventions are used: for a vector v in two-dimensional Minkowski space 
v'^Vp, = v+v- where v± = ± ui- The antisymmetric tensor is defined by = 1, r° = (t"/2 
with a°- being the standard Pauli-matrices satisfying a°'a'' = S"-'' + ie^'"^a'^, a,b,c— 1, 2, 3 and 
summation is implied over the repeated indices. 
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implying integrability (in the sense that an infinite set conserved quantities ex- 
ists). Next in Subsection 2.2 the axion model is demonstrated to be equivalent to 
the rather well studied SU(2)xU(l) symmetric anisotropic cr-model by a simple 
target-space duality transformation. This observation leads to a new Lax pair. 
Section 3 is devoted to perturbative investigations of the two models. First a brief 
renormalization group analysis of the anisotropic model is carried out, emphasiz- 
ing the important property of asymptotic freedom. Then the two loop /5-functions 
are compared in the two models with the conclusion that they are equivalent 
provided one takes into account a change of the renormalization scheme. In 
Subsections 3.2, and 3.3 the Legendre transformation of various (bosonic and 
fermionic type) finite density ground state energies (free energies) is computed 
in the anisotropic model. In Subsection 3.4 the main points of the analogous 
computation in the axion model are discussed. It is found that the free energies 
fully agree in perturbation theory. Since one compares physical quantities in the 
two (dual) models, no scheme dependence arises and therefore this provides good 
evidence for the validity of a quantum version of this T-duality. 
The anisotropic a-model is generally believed to be quantum integrable. Apply- 
ing the thermodynamical Bethe Ansatz one can also compute the free energies 
nonperturbatively from the bootstrap S'-matrix and the comparison with the 
perturbative results yields quite a stringent test of the proposed S'-matrix. This 
program is carried out in Section 4. 

First in Subsection 4.1, in the framework of the form factor bootstrap approach, 
the 2 — > 3 particle production amplitude in the 0(3) a-model with a constant 
^-term is shown to be non vanishing. To the best of my knowledge this is the first 
quantitative evidence for the (generally expected) quantum non-integrability of 
this model. To my mind this fact underlines the rather surprising nature of the 
quantum integrability of the axion model. Subsection 4.2 is devoted to the ther- 
modynamical Bethe Ansatz and in 4.3 the exact relation between the infrared 
and the ultraviolet mass scales (m/A-ratio) is given. Finally in Subsection 4.4 
the limit p —>■ (0(3) limit) is discussed, in particular it is explained how one 
recovers the known m/A-ratio of the 0(3) a-model. 

2 Lax pairs and T duality 

2.1 A Lax pair for the axion model 

The axion model belongs to a family of 0(3) symmetric classically integrable 
cr-models discovered in Ref. some of whose results I now recall. Introducing 
the matrix valued current 

= ^^^M-^'^ - ^^^^d-'n + l^d.n , (2.1) 
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where n = in°'a°', the equations of motion of ( p. . 1| ) can be written as: 

d^I^ = , d^I, - dj, = [I„ h] . (2.2) 



The standard form ( |2.2| ) of the equations of motion allows for the introduction 
of a Lax pair, 

U± = /± , (2.3) 

1 ± a; 

satisfying the zero curvature equation 

d,U,-d,U, = [U„U,], (2.4) 

for all values of the spectral parameter u. The current, J^, is closely related 
to the matrix valued Noether current, A/"^ = — ir^A/"", defined by 6C = d^e"'N^ 
corresponding to the symmetry transformation dn"" = e°'^'^e^n'^: 

I, = \U^ + e,,d-T, T=[^e-^)n. (2.5) 

Note that the trivially conserved part of {e^yd^T) is essential for the zero 
curvature equation (p^) to be satisfied. 



2.2 T-dual of the axion model 

Here I show that the axion model ( |1 . 1| ) is classically equivalent to a rather well 
studied a-model, the so-called anisotropic SU(2) principal a-model (PCM). The 
equivalence is a canonical one given by an Abelian target space (T-) duality trans- 
formation P, 1^). 

The Lagrangian of the anisotropic PCM can be written as: 

C^ = -^{LlL'^^ + gLlL^^], (2.6) 

where 

L, = G-^d,G = r''Ll, (2.7) 

and g is the parameter of anisotropy. The Lagrangian ( p.6|) can be interpreted 
as a deformation of the SU(2) xSU(2) (or 0(4)) symmetric nonlinear cr-model 
by the parameter g. For a generic value of g this (torsionless) model has an 
SU(2)lxU(1)r symmetry. At the classical level the anisotropic model, ( |2.6| ), 
interpolates between the SU(2)xSU(2) {g = 0) and 0(3) {g = —1) models. 

The equations of motion of the anisotropic principal model ( |2.6| ) can be writ- 
ten entirely in terms of the current L^: 

= , d^'Ll = -igL'^^Ll , d^'Ll = igL^^Ll . (2.8) 
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It is known that the equations of motion ( |2.8| ) admit a Lax representation 

i.e. there is a spectral parameter dependent current, V"^ = t"V^, satisfying the 

zero curvature equation (p.4|). This current can be written as: 



Vl'' = a±L'/, Vi = a±Ll, (2.9) 

where 



«± = ~1 2^ A ' ^± = ~1 2^ A ■ 2.10 

4 — gfo;"' ± 4a; 4 — guj'^ ± 4^^; 

To exhibit the classical T-duality transformation mapping the axion model 
( |1.1|) to the anisotropic PCM (|2.6| ) it is convenient to parametrize the n°-fields 
in Eq. (|1.1|) as 



An I 

n"*^ = sin -i? sin , = sin ■(9 cos , n^ = cos^9, 6 = — —Jl + gx-, 

^ (2.11) 

in terms of which the Lagrangian of the axion model ( |1.1| ) (after an integration 
by parts) becomes 

C = ^id^W^"^ + sin^ i^df^^d^"^ + d^xd^X + 2^Jl + g cos^e^'^df^xdu^} ■ (2.12) 
2A 

An Abelian T-duality transformation with respect to the S0(2)£i symmetry of 
the axion model ( p. . 1|) , corresponds to the following canonical transformation 0: 

/ ^ VI +^ / /o TON 

X = -TTTl^" = ' ^'-''^ 

where (and in the following) resp. Pa denote the canonical momenta conjugate 
to X resp. to its 'dual' a. In terms of the new variable, a, the dual Lagrangian 
turns out to be: 

£s = + {l + g cos^ §)d^,^d^'^ + (1 + ^) [d^ad'^a + 2 cos^d^ad^^] } . 

(2.14) 

A simple calculation shows that ( p.l4| ) is nothing but the Lagrangian ( |2.6|) , when 
parametrizing the SU(2) valued field, G, by the Euler angles 

Q = ^i^r^ ^Mr^ ^iar^ ^ (2.15) 

and taking into account the relations at the classical level between the couplings: 

A = A, ~g = g. (2.16) 

The observation that the axionic model is the T dual of Cy, also explains why 
9 should be interpreted as an angular variable. Indeed, as it has been shown in 
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Ref. [0 the Abehan T duahty transformation ( |2.13| ) maps the target space of the 
isotropic PCM, S^, the Abehan T duahty (|2.13|) into S"^ x S^. The arguments 



of Ref. [0] can be easily adopted to the present case with g > —1, and it is clear 
that in Eq. (|1 . 1| ) ra" parametrize an 5*^ and 6 parametrizes an S^. 

One can now use the classical T-duality transformation ( |2.13| ) to map the 
linear system of the axion model ( |2.2| ) to a new Lax pair for the anisotropic a- 
model (|2.6| ). It is given by Eq. ( p^) , where the current, J^, has to be replaced 
by 

= d^GG-' + g(Gr^G-')Ll - t^e^^d'' (Gt'G-') . (2.17) 



Eq. (|2.17|) is obtained from ( |2.1| ) by the T-duality transformation (p.l3|) . 

is related to the Noether current //^, corresponding to the manifest symmetry 

SG = —ie°'T°'G of (|2.6| ) and can be written analogously to A/'^^: 

4 = + ^f^'^d'f , f = -t^Gr'G-' . (2.18) 



It is clear that the new Lax pair ( p.3|) and the 'old' one, (|2.9|) , cannot be related 
by a gauge transformation since they have different pole structures as functions 
of the spectral variable, u. In the g —1 limit, the axion model reduces to 
the original 0(3) cr-model (decoupled from the 9 field), and the Lax pair ( p.3| ) 
becomes equivalent to that of Ref. 0, where it has been pointed out that the 
corresponding J^'s are ultralocal currents. I remark that the Lax pairs ( p.3|) and 



(2.9) correspond to (different) deformations of the usual Lax pairs of the principal 



chiral a-model, linear in d^GG ^ respectively G ^d^G. 



3 Quantum equivalence 
3.1 Perturbative investigations 

The canonical equivalence between the axion and the anisotropic PCM models 
naturally raises the question if there is also a quantum equivalence between them. 
Or putting it a slightly different way, one can ask if there is a quantum version 
of T-duality transformation ( |2.13| ). 

The equivalence of the models in perturbation theory (PT) is certainly a 
necessary condition for the validity of quantum T-duality. I start by recalling 
some renormalization properties of the anisotropic cr-model and the important 
property of asymptotic freedom. Then I present a quantitative comparison in PT 
between the two models by comparing their /5-functions. 

Due to its SU(2)lxU(1)r symmetry, the anisotropic cr-model ( pl6|) is renor- 
mahzable in perturbation theory. (Of course, in Eq. (|2^ ) (A, g) should be in- 
terpreted as the bare couplings (Acf^o)-) Using dimensional regularization, the 
relation between the bare parameters, (Ao,5'o), and the renormalized ones, (A, (7), 
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is written as: 



Xo = fi'XZx, and 1 + Qq = {1 + g) Zg . (3.1) 

In Eq. (|3.1|) the dimensional parameter, /x, is introduced as usual, to carry the 
mass dimension of Aq (which is dimensionless in 2D). The renormalization con- 
stants contain only simple poles in e = 2 — n whose residues can be calculated in 
perturbation theory. 

Physical quantities depend on the renormalized couplings A, g and the dimen- 
sional parameter, fi, in such a way that the action of the renormalization group 
(RG) operator 

d d d 

= ^ + /3a(A, g)g^ + /?,(A, g) - , (3.2) 

vanishes on them, where the /3-functions are obtained from the residues of the 
first order poles. The RG equations for the running couplings A, g can be written 
as 

^ = PxiXg), A(0) = A, (3.3) 

^ = Pg(X,g), m=9- 

One is interested in the 'large time' asymptotic behaviour of the couplings as 
functions of an external field (or chemical potential) h = /iq e*, with Hq fixed and 
t — > oo. In the anisotropic PCM the /3-functions are known to be of the form: 

(3a = -^{l -9 + Ap2(9) + ••■}, (3.4) 

where the two-loop /^-function coefficients are: 

1 - 2^ + 5g^ 1-g 
P2{g) = , and q-^[g) = — — . (3.5) 

Asymptotically free (AF) behaviour occurs when A — > together with g ^ gi 
as if: — i> oo, where gi is some constant. Analyzing Eqs. (|3.4| ), (|3.5|) , one finds three 
different AF solutions: 

1. ^1 = g = go = g = 0, (3.6) 

2. gi = -l g = go = g = -l, (3.7) 

3. gi = -l -l<g<0. (3.8) 



Solutions ( |3.6| ) resp. ( |3.7|) correspond to the 0(4) resp. the 0(3) a-model, while 
the solution (|3.8|) corresponds to the 'generic' case of the anisotropic PCM. 
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It is very convenient to introduce an exact RG invariant combination of the 
two couplings 

1 + git) , , 

p = 27rlim— 3.9 

and an effective /3-function for A(t) by 

/3eff(A,p)=/3A(A,r(A,p)), (3.10) 

expressing gif)^ in terms of the running couphng, A(t), and the RG invariant 
quantity p as git) = r(A(t) ,p). Using the perturbative result for r(A(t),p) 
one finds 

PAKp) = PAKp) = + + ■ ■ • • (3.11) 

Eq. ( |3.11|) allows for the introduction of an RG-invariant A-parameter in the MS 
scheme in the usual way: 



Ams = (A) (i^^) e^{i + oiX)} . (3.12) 

The effective coupling X^nih) is defined by the transcendental equation 

|^ + (f-l).„|^^.„^. (3,13) 

which is a function of the 'physical' parameter 

s = ln-^, (3.14) 

only. Moreover the running coupling can be expressed in terms of X^tfih) per- 
turbatively (in the sense that it is an infinite power series). The asymptotic 
expansion (for large s) of the effective coupling (containing terms oc In s) can be 
written as 

2n 7r(n-2), , , 

A.fl = — + M ng + --- . 3.15 

s 



After this detailed discussion of the renormahzation properties of the anisotropic 
PCM model I just briefly recall the main point for the axion model. As discussed 
in detail in Ref. the axion model is also renormalizable in the field theoretical 
sense and at one loop order the /3-functions of the couplings of the axion model, 

Pg, are simply obtained from /3a, Pg by the classical relation ( |2.16| ). At two 
loops, however, it has been found in [^] when using the background field method 
and dimensional regularization that the following perturbative redefinition of the 
couplings 

~X = X+^(l+g), ~g=g + A.{l+g)\ (3.16) 

(i.e. a change of scheme) is induced by the T-duality transformation. Taking into 
account Eqs. ( p.l6| ) the two loop /3-functions of the two models turn out to be 
equivalent. Alternatively, defining the renormahzation group invariant combina- 
tion of the two couplings analogously as in the PCM model Eq. (p.9|): 



p = 2vrlimiiM, (3.17) 
X(t) 

one finds 

p = p 

up to two loops . Thus as far as coupling constant renormahzation is concerned, 
the two models are equivalent, both are asymptotically free, and the actual value 
of p effects only the two loop coefficient. 



3.2 Perturbative tests of free energies 

Up to now quantum equivalence between the two dually related models has been 
tested by comparing simple renormahzation properties and the (two loop) /3- 
functions. The fact that the higher coefficients of the /3-functions are scheme 
dependent makes such a comparison more difficult and less conclusive. To lend 
some more credibility to the proposed perturbative quantum equivalence between 
the two models it is clearly desirable to compare physical quantities. For this 
purpose it is quite convenient to compute various zero temperature free ener- 
gies obtained by minimizing the Legendre transform of the Hamiltonian density 
coupled to some conserved currents 

'H = 'Ho-hiJi, H = J dxn = H -hiQi, (3.18) 

where hi are constant external fields having the interpretation of chemical po- 
tentials. In the limit of large fields {hi — >■ oo) corresponding to considering the 
system at large density, due to asymptotic freedom weak coupling perturbation 
theory (in A resp. A) is applicable. 
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The Legendre transform of the ground state energy density of the system must 
be of the form: 

J^{h) ^ 5f{h) ^ f{h) - /(O) = -h-'F.ih/k, Q) , (3.19) 

where A is a parameter of dimension mass. (Eq. ( |3.19D follows simply from di- 
mensional analysis.) In an asymptotically free theory the result of a perturbative 
computation of Fq (which is of course renormalization group invariant) is an 
asymptotic series in the running coupling, X{h/A) oc [ln(/i/A)]~^, where A is the 
usual renormalization group invariant combination of the (arbitrary) mass scale, 
fi, and the renormalized coupling, A. 

Since J-'{h) is a physical quantity, the results obtained in the two models 
should agree. Thus comparing them provides a further nontrivial check on the 
quantum equivalence between the axion and the anisotropic cr-model, hence also 
on the validity of the T-duality transformation at the quantum level. 

Since the axion field, 9, is actually an angle, its winding number (topological 
charge) can be non trivial. Therefore I present here the Legendre transforma- 
tion of the modified Hamiltonian (|3.18|) for a rather general case including also 
topological charges. 

Let us consider a general sigma model with torsion 

Co = \gABd^X%X'' + UABt^'d^X^d^X'', (3.20) 
and the following Ansatz for a set of conserved currents 

j; = cAX)d,x^ + e;:B\{x)d^x^ , (3.21) 

sufficiently general to include topological currents. The Legendre transformation 
of (|3.18|) yields the Lagrangian of the modified model which can be written as 

t = Co + h'Ji + hi'yC\C^^ . (3.22) 

In fact C can be obtained by gauging Cq i.e. by the substitution 

d^X"" ^ d,X^ + /i^5^oC^^ (3.23) 

when the antisymmetric field hAB is invariant (without compensating gauge trans- 
formation) under the symmetry transformation generated by the conserved cur- 
rents drii) . 

3.3 The anisotropic PCM model 

I recall next the main results of the perturbative computation of JF in the anisotropic 
PCM ( |2.6| ) first (for more details see Ref. [^). There are two conserved Noether 
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charges in the model ( p.6| ), Ql, Qr, corresponding to the U(1)lxU(1)r transfor- 
mation 

SG = ieicr^G - ienGa^ . (3.24) 

Introducing two chemical potentials coupled to the Ql resp. Qr, charges the 
Hamiltonian (|3.18|) takes the form: 

H = H^- hi^Qi^ - h^Qn . (3.25) 

Since the anisotropic PCM is torsionless the Legendre transformation of ( |3.25| ) 
can be simply found using the gauging procedure ( p.23|) . To actually calculate 
the ground state energy in PT one has to pass to the Euclidean field theory for- 
mulation in which case the chemical potentials hi can be interpreted as constant, 
imaginary gauge fields. The corresponding covariant derivative is then given 
by 

D2G = d2G + hLa'^G - huGa^ , D^G = d^G . (3.26) 

The gauged Lagrangian can be written as: 

C = Co + Ci + C2, (3.27) 

where Ci and £2 denote the terms linear and quadratic in the external fields, 
respectively. 

The computation of the free energy is based on the generating functional 



- / d"xf(h) 



lvGexp{- I (Tx {Co + C1 + C2)}. (3.28) 



In PT one expands Eq. ( 3.28|) in powers of the (bare) coupling Aq. I have also 
written the volume element as d^x, where n = 2 — e to indicate that dimensional 
regularization is employed in the perturbative calculations. 

In perturbation theory the action has to be expanded around a (stable) so- 
lution of the classical equations of motion. This amounts to find elements of 
SU(2), Gq, that correspond to (local) minima of the quadratic part of the gauged 
Lagrangian, 

r 2^ 2^(1 + ^?o) , ^hLhR{l + go)z 

L,2 — — —[^i+goz) \ , (_d.2yj 

Ao Ao Ao 

where 

^ = iTr {a^G- V^g} . (3.30) 
The first solution, referred to as 'bosonic', (BOS) is given as: 

Go = ^(1+^^x2), /i« = 0, (3.31) 
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and the other one considered here, which I call 'fermionic' (FER) is 



Go = ia\ (3.32) 

(The reason for the names 'bosonic' resp. 'fermionic' will become clear later on.) 

I start with the BOS case first. The leading, 0{\q^), term in perturbation 
theory is given by the potential energy, C2, at its minimum, which is immediately 
seen to be (zq = 0, = 0) 

^(-) = /:(-) = -^. (3.33) 

The 1-loop term leads to the following integral: 

^(0) _ /■ {'^ + 9o) pI - 9oP^ ,^ 

n J (27r)" (p2)2 _ ^g^hlp^ + 4(1 + g,)hlpl ' ^^"^ ^ 

For details concerning the actual computation, see Ref. 0, here I just quote the 
final RG improved perturbative result for the asymptotic expansion of the free 
energy for large values of the external fields: 

^Bos(/^) = + (1 - |) Ins + (ln2 - ^) + ■ ■ ■ I . (3.35) 

In the FER case the leading (classical) term of the free energy is 

^L'^) = C^-'^ = -^ii±M + hnf . (3.36) 

^0 

Note that JF(~^) for the diagonal charge depends only on the sum h = hi + hji. In 
the present case one needs to expand J-'{h) up to C(A^g.) in the effective coupling, 
which would at first sight necessitate a three-loop computation. In fact one 
obtains all the 0(Afg) terms from the one loop result alone! This 'mini miracle' 
is due to the (1 + g)"^ factor in front of the one- loop term together with the fact 
that (1 + g)"^ = 0(Afg.), implying that the one-loop term is already O(Afff). The 
result turns out to be: 



- "4^ ^« + 32^ 8^ 



Inp+ln 



2^^ 



+ 0(AL) . (3.37) 



Note the non-analytic contribution, InA^ff, in the last term of ( 3.37 ). The large s 
expansion in the FER case is finally given as 



1 , P 
Inp + - 

4 



+ ■ ■ ■ . (3.38) 



The above form of J-'{h) Eq. (p.37| ) explains the name 'fermionic' as its leading 
term corresponds to the free energy of p free fermions. The 0{l/s) terms can 
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be interpreted as interaction terms. In the 'bosonic' case the leading term in Eq. 
(|3.35|) is logarithmically divergent. 

Finally note that the one loop term also depends only on the sum hi + h^, 
just like the classical one (|3.36| ). This is not an accident. As shown in Ref. 0, 
in the FER case the free energy does depend only on the sum hi + to all 
orders in perturbation theory. In fact it is essential that J-'{hL, hpi) be a function 
oi h = h]^ + hfi in order to be able to match the perturbative result with the 
corresponding one computed by the non-perturbative Thermodynamical Bethe 
Ansatz for the DIAG case in Subsection 4.2. 



3.4 The axion model 

Next I outline the computation of the corresponding ground state energies to 
one loop order in the axion model ( |1.1| ), starting with the BOS case first. With 
the Euler angle parametrization of G ( p.l5|) the Ul(1) transformation, G i— >■ 



e**^"^ G of the anisotropic a-model ( |2.6D acts as a simple shift, (p{x) i— > (p{x) + 
K. The corresponding Noether charge, Ql, and its image under the T-duality 
transformation, Ql? are simply 

dxp^ , Q^ = J dxp^ , where p^ = ' ^ ^ ' 



since the canonical transformation implementing the T-duality mapping ( |2.13| ) 
effects only p^, x'; ct' and p^, leaving the other fields, y), ^9, p^, p^, unchanged. 

Since in the BOS case the h^B field in Eq. (|1 . 1| ) is invariant, one can simply 
'gauge' the Lagrangian of the axion model in an external {hi) field (see Eq. (|3.23|) . 
The classical ground state is found to be = x = 0, = 7r/2. (The corresponding 
solution of the anisotropic a-model is given by = a = 0, = vr/2.) 

Expanding the (Euclidean) Lagrangian (after suitable rescalings, etc.) one 
obtains 

_ 2/;2 1 

C= -^ + -mMm^ + o{\), (3.39) 
Ao 2 

-d''+Ahl 2/iLv^r+^e^25^ 

M=\ -d^ I , (3.40) 



where 



|2 



and m = (t?, 99, x)- (Ao, ^0 denote the bare coupling and parameter of the axion 
model). In Eq. (|3.40|) the e tensor has been explicitly kept, as it requires a careful 
definition in n = 2 — e dimensions used to regularize the momentum integrals. 



It is convenient to adopt the definition of Ref. ||T0|, where this antisymmetric 
tensor corresponds to an almost complex structure: e^y = —e^^, e^^e^^j = dy^j. 
The one loop quantum corrections to the classical ground state (the first term in 
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Eq. ( p.39| )) require the calculation of a functional determinant, leading to 

= / fn -4 I AU2f-2 ' Pt^ = ^t^'^P'^ ' 3.41 

n J {2tc)'' p'^ + Ahi{pi - gopi) 



To evaluate (|3.41| ) one should apply the modified dimensional regularization of 
Ref. as p2 = ^2uPu plays a distinguished role and it is kept as a one dimensional 
variable. In fact for our purposes it is sufficient to calculate the difference J-'{h) — 
^j:{h), where J^j]{h) is the corresponding determinant in the anisotropic a-model 
(Eq. (3.12) in Ref. 0). Since both J-'{h) and J-'-s^h) are already the first quantum 
corrections to the classical expressions one may set g = g (and make no distinction 
between bare and renormalized g's) when computing their difference to lowest 
order and one ends up with 

Hh)-Mh) = -^(^+9)J J^^^J^ = -^a + 9M9), (3.42) 

where Ni = q"^ + qf — gq2, N2 = q'^ + q2 — gqf- Although the integrand yielding 
w{g) is antisymmetric under qi <-> q2, the integral is divergent by power counting 
for n = 2, i.e. it must be computed in n = 2 — e dimensions. Its derivative, 
w'{g), is, however, convergent by power counting and it has also an antisymmetric 
integrand, therefore this latter may be evaluated in = 2 dimensions giving 
w'{g) = 0. Then to compute w{g) one may choose e.g. the point g = —1: 

r d'^q ql-ql _ n-l-l r d'^q q^ _ 1 
' J (27r)" (g2 + 1)2 n J (27r)" (g^ + 1)2 4vr' ^ ' ' 

where writing the second equality, it has been used that qi is an — 1-dimensional 
variable, while q2 is 1-dimensional. From ( |3.43[ ) one finds that after taking into 



account the change of the renormalization scheme ( p.l6 ), in PT the free energy 



densities of the two models ( |1.1[ ) and ( |2.6| ) do indeed coincide for the BOS case. 



I omit the details of the calculations in the PER case (see Ref. [|TT] for de- 
tails), suffice it to say that this free energy is also in perfect agreement with the 
corresponding result in the anisotropic PCM model. 

I would like to point out here that the canonical transformation connecting C^, 
and £, maps Qr to a purely topological charge Qr of the axion model, which is 
completely different from the Noether charge of the 'manifest' SO(2)0 symmetry 
of the Lagrangian ( |1 . Ij ) . To illustrate this I quote the value of the classical free 
energy density corresponding to the Noether charge of the SO (2)51 symmetry: 
-^Imin = ~2/i|/A, very different indeed from Eq. ( p. 361) with /il = 0. 



In conclusion the free energy densities fully agree for the BOS and PER cases 
in both models up to first order in RG improved PT, providing a rather convincing 
evidence for the validity of T-duality at the quantum level. 
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4 Quantum integrability 



It is by now generally accepted that the anisotropic PCM model (|2.6| ) is integrable 
at the quantum level ^ in the sense that there is no particle production. 

Assuming that the spectrum contains two massive doublets (kinks) their scatter- 
ing is described by the tensor product of an SU(2)xU(l) symmetric solution of 
the bootstrap S'-matrix equations: 

= S(°°)(^)®5(p)(^), (4.1) 

where S^p\9) denotes the Sine-Gordon (SG) S'-matrix. Let me point out here 
that denoting the parameter in the 5* matrix in the same way as the RG invariant 
quantity in Eq. ( p.9| ),(p), is intentional. As will become clear from the results 
later on, one should actually identify them. (In fact this identification is seen to 
be consistent at least up to two loops in PT.) 

In the limit p —>■ oo (corresponding to g = —1) the symmetry of the S'-matrix 
([4.1|) increases to SU(2)xSU(2). For the range of the parameter < p < 1 
in addition to the kinks bound states (breathers) also appear in the spectrum, 
transforming as 3 + 1 under SU(2). 

In this section I shall explore some consequences of assuming the validity of 
the duality transformation at the full quantum level between the two theories. 
This immediately implies the absence of particle production in the axion model 
(|1.1|) too, and that its two particle S'-matrix is given by Eq. ( |4.1| ). Note that the 
global symmetries of the axion model, SO(3)xSO(2)e, cannot be identified with 
those of the S matrix ( |4.1| ). As it will be convincingly shown the U(l) symmetry 
of the S matrix ( [4.1|) corresponds to a hidden 'topological' U(l) of the axion 
model, the dual of the U(1)r Noether symmetry of the anisotropic PCM under 
the T-duality transformation. 

The following (rather heuristic) consideration might be useful to give some 
insight into the connection between the SG and the axion model. Integrating 
out the 0(3) fields, n"", in some generating functional of the theory ( |1 . 1| ) one 
obtains a non-vanishing effective potential for the 6 field. Since 6 is 27r-periodic 
the effective potential must be also periodic. The effective theory of the 6 field 
is therefore expected to show some features similar to the Sine-Gordon model, 
with a periodic potential and corresponding topological current = e^j^^d'^ 6/271. 
(Note that 6 being an angular variable makes it difficult to integrate it out in the 
functional integral in spite of C being only quadratic in 6.) 

4.1 Particle production in the 0(3) model with a ^-term 

At first sight the proposed quantum integrability of the axion model seems to 
be rather questionable at least, as it is generally believed that the 0(3) cr-model 
with a constant 6'-term is not quantum integrable, except for the special value 
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6 = n (despite the fact that the 6'-term, being a total derivative, does not 
change the classical physics of the model). 

In the framework of the form-factor bootstrap approach one can actually 
show that the 6 term mediates particle production, indeed [|lT] . To illustrate this 



important fact, let us express the 2 — > 3 particle production amplitude to lowest 
order in 6 as: 

{p, b; p', b'; /, b"\q, a; q', a')^e) = {27rf 1 9 6'^'\p + p' + p" - q - q') 

■ {p, b; p', b'; p", b"\T{0)\q, a; q', a')(o) + 0(6') , (4.2) 

where in the first line the amplitude is in the 0(3) model with a 6'-term, while in 
the second line the matrix element of the topological charge density operator T 
is to be calculated in the original 0(3) cr-model (with ^ = 0). In other words one 
computes the 2 — > 3 particle production amplitude to first order in 6. 

For the sake of simplicity let us consider the following special kinematical 
configuration: the incoming particles have momenta qi = Q and q'^ = —Q, 
whereas the produced (outgoing) three particles have momenta Pi = Q', p[ = 
and p'l = —Q' respectively. Here Q' can easily be expressed in terms of Q and 
the kink mass M using energy conservation. For large Q, using the results of Ref. 



ril, one finds 



{p,b;p\b';p'',b''\Tmq,a;q\a')^,)^nl^^^j^ . (4.3) 



Eq. ([4.3|) shows that already to first order in 9, the 2^3 particle production 
amplitude is different from zero. Thus at least for small values of 9, the introduc- 
tion of this term destroys quantum integrability of the 0(3) cx-model, indeed. To 
the best of my knowledge this is the first quantitative evidence for the quantum 
non- integrability of the 0(3) cr-model with a constant 6'-term. 



4.2 Free energy from the Thermodynamical Bethe Ansatz 

The free energy of a system of particles of mass m is of the form: 

J^{h) = -h^Fo{h/m,Q). (4.4) 

The Thermodynamical Bethe Ansatz (TBA) method based on the S'-matrix of 
the particles leads to a set of integral equations for the free energy whose solution 
can be found in many cases as an asymptotic series expansion in /i/m ^ 1. Since 
^ is a physical quantity, the results obtained by TBA method and by perturbation 
theory should agree. Thus by comparing the asymptotic series of FQ{h/m,Q) 
in h/m with the one of FQ{h/A,Q) in h/A obtained in PT provides a rather 
stringent consistency check on the S-matrix (and also on the self-consistency of 
the hypothesis used in the course of the calculation). Moreover one obtains the 
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exact m/A ratio, a rather important non-perturbative parameter of the theory 
which can be measured e.g. by lattice simulations. 

For definiteness we choose hi^,hii> and assume p > 1, implying that there 
are two doublets in the spectrum. The charges of the particles are normalized as: 

(1,1), (-1,1), (1,-1) (-1,-1). (4.5) 

One can distinguish between the three possible types of finite density ground 
states depending on /il, /?-r as: 

1. /iL , /iR > , (DIAG) 

2. = , /iR = /i > , (RIGHT) 

3. = /i > , /iR = . (LEFT) 

• The DIAG case. 

As the system ( |3.25|) is considered with hi^/m, hn/m ^ 1, it is clear that par- 
ticles of charge (1,1) condense into the vacuum. It is less clear what other kind 
of particles (necessarily with a smaller charge/mass ratio) will appear in the vac- 
uum state. I shall assume that the vacuum consist of only particles of charge 
(1,1). This seemingly radical assumption has apparently worked in all analogous 
examples studied so far, and it greatly simplifies the solution of TEA equations. 

The scattering phase of the (1,1) type particle, S{9), can be easily found from 
the S-matrix: 

6(6) = 6^(9) + 6pie) . (4.6) 

When the ground state is assumed to contain only particles of charge (1,1) the 
calculation of the free energy using the TEA method reduces to the solution of a 
single integral equation ( [4.8|) whose kernel is given by the logarithmic derivative 
of the relevant S'-matrix element. 

The integral equation in this case is given as follows: 

e(e) - r de'Kie - e')de') = h-m cosh e , (4.8) 

J-B 

where h = hi^ + hn, together with the boundary condition e(±i?) = 0. In terms 
of the solution of Eq. ( |4.8| ) the free energy is 

Sf(h) = -— d9coshi9)e{0). (4.9) 

ZTT J-B 

The asymptotic series for the free energy up to 0(lnt/t^) terms reads as: 



6Uh) = --p 

71 



2t 4 t2 + ^2 + ^3 ' 



(4.10) 
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where t = \n{h/m) and Ad is a constant which determines the m/A ratio: 



2 



lnr(l + |) + - ln2) - 1 + 31n2 - InTr 



(4.11) 



Let me remark here that from the classical term in Eq. ( |4.1CI| ), — /i^p/vr, one can 
immediately read off the level of an underlying ultraviolet (UV) current algebra, 
k, as k = p [ p!5|] . 

• The RIGHT case. 

In analogy to the assumption made in the diagonal case one would now expect 
the ground state to consist of a mixture of particles of the same charge with 
respect to U(1)r, i.e. those of charge (1,1) and (-1,1) with equal densities. In 
Ref. [|I^ the result of the diagonalisation of the pertinent coupled TEA system 



has been given and one finds (somewhat surprisingly) that 

Sfnih) = 6Uh) , (4.12) 

where h = /ir. So the free energies are precisely the same as functions of the 
effective chemical potentials in these two apparently rather different cases. 

• The LEFT case. 

Analogously to the RIGHT case the ground state is expected to consist of an 
equal density mixture of particles of the same charge, this time with respect to 
U(1)l, i.e. those of charge (1,1) and (1,-1). Here I just quote the asymptotic 
series of the free energy density: 

^/l(^) = - V + (1 - |) Int + + . . .) , (4.13) 
where the constant Al, determining the ratio m/A is: 

^L=(3-|)ln2 + ^-lnvr + lnr(|) . (4.14) 



• The LEFT case for p < 1. 

In the previous subsections it has been assumed that the parameter p > 1. In 
order to clarify what happens when p becomes smaller than one, I recall that then 
bound states (breathers) appear in the spectrum (for p < 1 one enters simply the 
attractive regime). The SG mass spectrum is given by the well known formula: 

TTpr 1 
rrir = 2m sin , r = 1, 2 . . . < - , (4-15) 

where m denotes the mass of the SG kinks. (Note that m also gives the mass of 
the SU(2) doublets.) The appearance of these new particles implies that as soon 
as p becomes smaller than 1 the vacuum changes as then it becomes energetically 
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favorable for the r = 1 charge (2,0) breather to condense. This charge (2,0) 
breather can be interpreted as the bound state of the (1,1) and (1,-1) particles 
(kinks). In fact as this particle has the highest charge/mass ratio, one can again 
assume (as in the DIAG case) that the true vacuum consists only of the conden- 
sate of the r = 1 charge (2,0) breather. Calculating the breather-breather phase 



shift by the bootstrap-fusion method [|T6| and solving the integral equation ( |4.8| ) 



yields the free energy density of the breather condensate, 5fB{h-B,mi), where 
mi = 2m sin ^ denotes the breather mass. One finds that 5/b(/'-B;^i) becomes 
identical to the free energy of the charge (1, ±1) kink mixture, Sfi^{h, m) analyti- 
cally continued for p < 1. This is very reassuring as the result for Sfi^{h, A) from 
perturbation theory turns out to be continuous for the full range < p < oo, 
which is certainly to be expected. 

4.3 The m/Ay^ ratio 

Having calculated the asymptotic form of the free energy both in perturbation 
theory and with the TBA method, by comparing them (with t = s — ln(m/A(^)) 
one obtains the relation between the mass of the doublet particles, m, and A^. 
In both the BOS (LEFT) resp. FER (RIGHT,DIAG) cases the comparison of 
Eqs. i fn^ and QO^ ) resp. and {^^^ leads to the result 



- _o3-.,l-ir(i + p/2)_ ^^^^g^ 



I would like to emphasize that it is already a very nontrivial check on the overall 
consistency of the assumptions that all of the expansion coefficients, for both the 
BOS and the FER cases, agree. 

To illustrate how nontrivial this consistency is, let me point out the following 
(apparent) paradox. In Eq. ( [4.10D the ratio of the coefficients of the Int/t^ term, 
— and the —1/t term, p/2, is = — p/2, while in Eq. ( [4.13| ) the ratio of the 



coefficients of the Int term, 1 — p/2, and that of the t term is Tl = 1 — p/2. This 
is quite remarkable, since on very general grounds in a theory with a coupling. A, 
the perturbative result for a fermionic type free energy is expected to be of the 
form 

5/p.(/.) = -^ (l-h^ + OiA , (4.17) 
while in the same theory for a bosonic case it is expected to be given by 

SfBosih) = (^j^ + const. + Oixij , (4.18) 

where fo, /i, bo are constants and X{h) is the running coupling. Since from the 
renormalization group 

^ /3oln^ + ^lnln^ + 0(-i-), (4.19) 



X{h) X (3o A 'VA' 
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one expects rper = ''"Bos = Pi/Po- As in the present case 7^ Tl, it seems impos- 
sible to match Eq. ( [4 .101) and Eq. ( [4.13|) simultaneously with the corresponding 
perturbative expansion. The resolution of this paradox lies in the presence of the 
non-analytic contribution, In Aeff, in the last term of Eq. (|3.37|) . This non-analytic 
term in Aeff explains why the coefficient of the In s term differs for the BOS and 
the PER cases. Such non-analytic terms can only occur in perturbation theory 
in the presence of several couplings. 

A further, very stringent consistency check is that the m/A ratio obtained in 
the BOS case is exactly the same as the one obtained in the FER case. Finally 
by comparing the corresponding free energies of the breather condensate one 
obtains the m/A ratio in Eq. ([4.16|) . This is clearly an additional check on our 
result and on the mutual consistency of the hypothesises made in the course of 
the calculation. 



4.4 The p ^ limit 

The p — > limit (i.e. the running deformation parameter, g{h) —1), is expected 
to give the 0(3) a-model (and a decoupled free field) [0. In fact from Eq. (|4.15| ) 
one sees that m/mi — >■ cx) as p — i.e. the kinks 'disappear' from the spectrum. 
Furthermore the higher breathers with r > 1 also 'disappear' from the spectrum, 
but for a different reason. The binding energy of the r > 1 breathers (which 
can be considered as bound states of the r = 1 ones) tends to zero giving rise to 
zero energy bound states. Therefore the limiting spectrum consists of 4 massive 
particles, transforming as (3l + 1l, 1r) under SU(2)lxU(1)r. Our results lend 
some additional support to this expectation, as the mi/Aj^ = 8/e ratio in the 
0(3) cr-model calculated previously (where now mi stays for the mass of the 
triplet in the 0(3) model) is beautifully reproduced in the p — > limit. This 
can be immediately seen from Eq. ( 4.16 ) using the fact that in the p —>■ limit 
nil ~^ npm. 



5 Conclusions 

In conclusion one can say that there is complete consistency between the TBA 
and the perturbative calculations in both models, providing good evidence for 
the validity of the proposed S matrix (|4.1| ) for them. Since the effective cou- 
pling (|3.11|) is identical in the two models the m/Ajj^ ratio Eq. ( [4.16| ) stays also 
the same. It might be worth emphasizing once again that these results provide 
also nonperturbative evidence for the validity of quantum T-duality. In particu- 
lar one has some quantitative evidence for the identification between the U(1)r 
symmetry of the proposed S matrix and the topological current of the axion 
model (corresponding to the canonical image of the U(1)r Noether current of the 
anisotropic PCM and being quite distinct from the manifest S0(2)e symmetry). 
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The above findings provide good motivation to study both models by lat- 
tice Monte-Carlo simulations. Measuring the m/ ratios (A^, denoting the A- 
parameter of the lattice regularized theory) would provide us with a completely 
non-perturbative way of testing quantum T-duality. 
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